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Abstract 

Magnetic resonant coupling (MRC) is an efficient method for realizing the near-field wireless power 
transfer (WPT). The use of multiple transmitters (TXs) each with one coil can be applied to enhance the 
WPT performance by focusing the magnetic fields from all TX coils in a beam toward the receiver (RX) 
coil, termed as “magnetic beamforming”. In this paper, we study the optimal magnetic beamforming for an 
MRC-WPT system with multiple TXs and a single RX. We formulate an optimization problem to jointly 
design the currents fiowing through different TXs so as to minimize the total power drawn from their voltage 
sources, subject to the minimum power required by the RX load as well as the TXs’ constraints on the peak 
voltage and current. For the special case of identical TX resistances and neglecting all TXs’ constraints 
on the peak voltage and current, we show that the optimal current magnitude of each TX is proportional 
to the mutual inductance between its TX coil and the RX coil. In general, the problem is a non-convex 
quadratically constrained quadratic programming (QCQP) problem, which is reformulated as a semidefinite 
programming (SDP) problem. We show that its semidefinite relaxation (SDR) is tight. Numerical results 
show that magnetic beamforming significantly enhances the deliverable power as well as the WPT efficiency 
over the uncoordinated benchmark scheme of equal current allocation. 

Index Terms 

Wireless power transfer, magnetic resonant coupling, magnetic beamforming, semidefinite relaxation 

I. Introduction 

Near-filed wireless power transfer (WPT) has been attraeting growing interest, due to its high 
effleieney for power transmission. Near-held WPT is realized by induetive eoupling (IC) ifTl. If2l for 
short-range applieations in eentimeters, or magnetie resonant eoupling (MRC) S, [!4j| for mid-range 
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applications up to a couple of meters. With the short-range WPT teehnology reaehing the stage of 
eommereial use, the mid-range WPT teehnology has been gathering momentum in the last deeade. 
The reeent progress on mid-range WPT is referred to the review paper Q and the referenees therein. 

The MRC-WPT system with multiple transmitters (TXs) and/or multiple reeeivers (RXs) has been 
studied in the literature [|5l|-[l8l. The MRC-WPT system with only two TXs and one RX is studied 
in [|5||, im, while the obtained results eannot be readily extended to the ease of more than two TXs. 
Reeently, an “Magnetie MIMO” eharging system [|3 ean eharge a phone inside of a user’s poeket 
40em away from the array of TX eoils, independently of the phone’s orientation. For an MRC-WPT 
system with multiple RXs, the load resistanees of the RXs are jointly optimized in [|8l to minimize 
the total transmit power and address the “near-far” fairness problem. Deploying multiple TXs ean 
help foeus the magnetie fields on the RX [| 7 ]|, in a manner analogous to beamforming in far-field 
wireless eommunieations Q. However, to our best knowledge, there has been no prior work that 
designs the magnetie beamforming from a signal proeessing and optimization perspeetive, for an 
MRC-WPT system with arbitrary number of TXs, whieh thus motivates our work. 

In this paper, as shown in Fig. [U we eonsider an MRC-WPT system with a single RX and 
multiple TXs of whieh their eurrents (or equivalently souree voltages) ean be adjusted sueh that 
the magnetie fields are eonstruetively eombined at the RX, thus aehieving a magnetie beamforming 
gain. We formulate a problem to minimize the total power drawn from the voltage sourees of all TXs 
by designing the eurrents flowing through different TXs, subjeet to the minimum power required 
by the RX load as well as the TXs’ eonstraints on the peak voltage and eurrent. For the speeial 
ease of identieal TX resistanees and negleeting the TXs’ eonstraints, the optimal eurrent magnitude 
of eaeh TX is shown to be proportional to the mutual induetanee between its TX eoil and the 
RX eoil. In general, our formulated problem is a non-eonvex quadratieally eonstrained quadratie 
programming (QCQP) problem. It is reeast as a semidefinite programming (SDP) problem. We show 
that its semidefinite relaxation (SDR) [fT^ . ifTTll is tight, i.e., the existing optimal solution to the 
SDR problem is always rank-one. The optimal solution to the QCQP problem ean thus be obtained 
via standard eonvex optimization methods IfT^ . Numerieal results show that magnetie beamforming 
signifieantly enhanees the deliverable power and the WPT effieieney over the benehmark seheme 
of equal eurrent alloeation. 
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II. System Model 

As shown in Fig. [H we consider an MRC-WPT system with > 1 TXs each equipped with a 
single coil, indexed by n, and a single RX with a coil, index as 0. Each TX n is connected to a 
stable energy source supplying sinusoidal voltage over time given by Vn{t) = with the 

complex Vn denoting an adjustable voltage and ze > 0 denoting its operating angular frequency. Let 
in{t) = denote the steady-state current flowing through TX n, with the complex current 

in- This current produces a time-varying magnetic flux in the n-th TX coil, which passes through the 
RX coil and induces time-varying currents in it. Let io(^) = be the steady-state current 

at the RX, with the complex current zq. 

We use Mno and to denote the mutual inductance between the n-th TX coil and the RX 

coil, and the mutual inductance between the rzi-th TX coil and the n 2 -th TX coil, respectively. Each 
mutual inductance plays the role of magnetic channel between a pair of coils, and is a positive 
or negative real number depending on the physical characteristics, the relative distance, and the 
orientations of a pair of coils, etc. (3. 



Fig. 1: System model of an MRC-WPT system 

We denote the self-inductance and the capacity of the n-th TX coil (RX coil) by > 0 (Lq > 0) 
and Cn > 0 (Cq > 0), respectively. The capacitors are chosen such that all TXs and the RX have 
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the same resonant angular frequency w. We use > 0 to denote the sum of the source resistance 
and the internal parasitic resistance of the n-th TX. The resistance of the RX, denoted by tq > 0, 
consists of the parasitic resistance rp_o > 0 and the load resistance r| q > 0, i.e., tq = rp o + r| q. The 
load is assumed to be purely resistive. 

In this paper, we assume that all electrical parameters are fixed and known by the central controller. 
For analytical convenience, we treat the complex currents z„’s flowing through TXs as design 
parameteriil which are adjusted to realize magnetic beamforming. 

By applying Kirchhoff’s circuit law to the RX, we obtain its current zq as follows 


N 



( 1 ) 


Define the vector of mutual inductances between the RX coil and all TX coils as m = [Mqi M 02 • • • 
Mqn]^. From O, the power delivered to the RX load is 



• M M • 

1 mm 1 . 


( 2 ) 


By applying Kirchhoff’s circuit law to each TX n, we obtain its source voltage as follows 

N 

Vn = rJn + jw ^ Mnkik “ jwMnOlo 
k=l,^n 



( 3 ) 


We define the iV-order square matrix B as follows 


B = B + jB 


( 4 ) 


where the elements of the real-part matrix B and the imaginary-part matrix B are given by 



( 5 ) 


'in practice, it is more convenient to use a voltage source rather than a current source. However, by using standard circuit theory, 
one can easily compute the source voltages Vn’s that generate the the required currents i„’s. 
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Bnk 


0 , 

-wMnk 5 


if k = n 
otherwise 


( 6 ) 


The matriees B, B and B are symmetrie, due to the faet that Mnk = Mkm Vn ^ k. Denote the n-th 
eolumn of the matriees B, B, B by b„, b„, b„, respeetively. Moreover, the matrix B is positive 
semidefinite (PSD), as it ean be rewritten as follows 


B = R 




To 


(7) 


The souree voltage of eaeh TX n ean be equivalently rewritten as 


Vn = b^i. 


( 8 ) 


From dH) and ([8]), the total power drawn from the sourees of all TXs, denoted by p, is derived as 
follows 




-i^Bi. 

2 


(9) 


Remark 1. From (|5]) and dH), we observe that the total power drawn from all TXs’ voltage sourees 
depends on the mutual induetanee Mno between the eoils of eaeh TX n and the RX, but it is 
independent of the mutual induetanee M^k between any pair of TX eoils. 


III. Problem Formulation 

In this seetion, we formulate a problem to minimize the total power drawn from the voltage 
sourees of all TXs by jointly designing the eurrents i flowing through TXs, subjeet to praetieal 
eonstraints of the MRC-WPT system. Partieularly, we eonsider the following eonstraints: first, the 
power delivered to the RX load should exeeed a given power level /3o > 0, i.e., po > /So; seeond, the 
maximally allowable amplitude of the souree voltage Vn is I 4 . i.e., |u„| < 14 ; third, the maximally 
allowable amplitude of the eurrent 4 is An, i.e., |4| < An. Let Q„ be the rank-one matrix with the 
n-th diagonal element as one and all other elements as zero. From d2l), dS]) and dl]), the problem is 
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thus formulated as follows 


(PI): 

min -i^Bi 

iec^ 2 


(10a) 

s. t. 

2^2 1 mm i>/^o 


(10b) 


i^b„b^i<V;^ n = 

1, 2, ••• N 

(10c) 


i^Qni < Al, n = 1, 

2, ■■■ N 

(lOd) 


(PI) is a complex-valued non-convex QCQP problem ITT^ . Although solving such a problem is 
nontrivial in general ifTOll . ifTTI . we obtain the optimal solution to (PI) in Section HVl 


IV. Optimal Solutions 

A. Optimal Solution to (PI) without Constraints (IIPel) and (llOdI) 

In this subsection, we consider the simplified (PI) only with delivered load power constraint (IIObi) 
but without voltage and current constraints given in (IIPel) and (llOdI) . respectively, to get useful 
insights on magnetic beamforming. We observe that the real-part currents i and the imaginary-part 
currents i contribute in the same way to the total TX power in (II Pal) and the delivered load power 

in (IIObi) , since both B and mm^ are symmetric matrices. As a result, we can set i = 0 and design 

i only, i.e., we need to solve 

(P2) : min -pBi (11a) 

s, t, ( 11 b) 

2i-5 

The optimal solution to (P2) is given as follows. 


Theorem 1. The optimal solution to (P2) is i* = aui, where a is a constant such that the constraint 
(II Ibl) holds with equality, and Ui is the eigenvector associated with the minimum eigenvalue, denoted 
by 7 i, of the matrix 


T 


R 


w 


’{ro 


P,o H 

-mm , 


( 12 ) 


where v* is chosen such that 71 = 0 . 
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Specifically, for the case of identical TX resistances, i.e., R = rl, the optimal currents to (F2) is 


F = 


am 


™ 2 


Proof: We construct the Lagrange function of (P2) as 


L{i,v) = ^i^Bi + n ( /3o 
Then, the Lagrange dual function is given by 


-75—1 mm 1 

2rl 


L{v) = fiov + inf ( ^i^Bi - — 

te^mm^ 

n.o + Tpfl 


i“mm^i 


= fiov+ini 


P ( R + 




2rl 


TT— 

-i‘mm 1 


1 - / '■■'2 
= fiov+ini -Pf R4 


U)2((l-'i;)ri,o + rp,o) h , 

-- - —mm 1. 


(13) 


(14) 


(15) 


To obtain the best lower bound on the optimal objective value of (P2), the dual variable v should be 
optimized to maximize the Lagrange dual (fTSl) . For dual feasibility, the Lagrange dual function (fTSl) 
should be bounded below. For convenience, we denote the following singular-value-decomposition 
(SVD) 


R + 


W 2 ((l 


v)r\.o + rp,o) H 

—FT--—mm 


uru^, 


(16) 


where the matrix U = [ui U2 ■ ■ ■ u^v] is orthogonal, and P = diagjyi, • • • ,7Ar}, with 71 < 72 < 
• • • < 7 Af- For the case of arbitrary transmitter resistances, the Lagrangian in (fT4l) is bounded below 
in i and the Lagrange dual function (fTSl) is maximized, only when v is chosen as v* such that 

7i = 0. 

Moreover, we observe that the objective (lllal) is minimized when the constraint (IIIbl) holds with 
equality, since both B and the matrix mm^ are PSD. Hence, the optimal current can be written as 
i* = aui, where ui is the eigenvector associated with the eigenvalue 71 = 0 , and a is a constant 
such that the constraint (lllbl) holds with equality. 

For the case of identical transmitter resistance, i.e., R = rijv, from the isometric property of the 
identity matrix I^r, the diagonal matrix P is given by 


P = diag 
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and the eigenvector ui = and u„, Vn > 2, are arbitrarily orthogonal vectors constructed 


by methods such as Gram—Schmidt method. It is standard to show that the Lagrangian in (fT4l) is 
bounded below in i and the Lagrange dual function (fTSl) is maximized, only when v is chosen such 
that the first eigenvalue is zero, i.e., the optimal dual variable is 


n* = 1 + 


rr. 


r\nW^ 


0 


(17) 


and the optimal current is given by 


am 
1 = r— 

where a is a constant such that (II Ibl) holds with equality. 

In fact, all TX currents can be rotated by arbitrarily common phase. 


(18) 


Remark 2. Theorem [U implies that for the case of identical TX resistances, the optimal current 
magnitude of each TX n is proportional to the mutual inductance Mo„ between the RX and TX n. 
This is analogous to the traditional maximum-ratio-transmission (MRT) beamforming in wireless 
communications [Q. The magnetic beamforming also differs from the traditional beamforming in 
wireless communications. The former operates over the near-field magnetic flux, and only the TX 
current magnitudes are adjusted according to the real magnetic channels (i.e., positive or negative 
mutual inductances) [| 71 |; while the latter operates over the far-field propagating waves, and the 
amplitude as well as the phase of the signal at each TX are adjusted according to complex wireless 
channels including amplitude fluctuation and phase shift due to prorogation delay lfT3l . 


B. Optimal Solution to (PI) 

Define X = ii^, M = mm^, and B„ = b„b^. Thus, (PI) can be equivalently rewritten as the 
following SDP problem 


(Pl-SDP) : min - Tr (BX) (19a) 

s. t. Tr (MX) > (19b) 

Tr(B„X)<K^ n = l, ••• iV (19c) 

Tr(Q„X)<A2, n = l, iV (19d) 

X 0, rank (X) = 1 (19e) 

September 24, 2015 DRAFT 








where X ^ 0 indieates that X is PSD. 

As well known, the rank eonstraint in (I19el) is non-eonvex ifTOl . [fTTI . By ignoring the rank-one 
eonstraint in (I19el) . we obtain the SDR of (PI—SDP), denoted by (PI—SDR), whieh is eonvex. 
Moreover, we have the following theorem on (PI—SDR) and the optimal solution to (PI). 

Theorem 2. The SDR of (PI—SDP) is tight, i.e., the existing optimal solution X* to (PI—SDR) 
is always rank-one (i.e., X* = i* (i*)'^). The optimal solution to (PI) is i*. 

Proof: Let A > 0, p = (pi, • • • , pat)^, Vp„ > 0, p = (pi, • • • , piv), Vp„ > 0, be the dual 
variables corresponding to the constraint(s) given in (IIObi) . (llOcI) . and (llOdI) . respectively. Let the 
matrix S 0 be the dual variable corresponding to the constraint X 0 in (I19el) . The Lagrangian 
of (PI—SDR) is then written as 

L(X,A,p,S) = ilV (BX) - A (t1'(MX) - + 

2 V ) 

N N 

J^p^Tr (B„X)-+5^p4Tr (Q„X) _Tr (SX). (20) 

n=l n=l 

Let X*, A*, p*, p*, and S* be the optimal primal variables and dual variables, respectively. More¬ 
over, the Karush—Kuhn—Tucker (KKT) conditions are given by 

.. N N 

VxP(X*, A*, pf S*) = -B-yM + J2 ^^nQn - S* = 0. (21) 

n=l n=l 

S*X* = 0. (22) 


Next, by postmultiplexing (|2TI) by X* and substituting (l22l) into the resulting equation, we obtain 

N N 


1 


-BX* - A*MX* + p:B„X* + p^QnX* = 0. 


(23) 


n=l 


n=l 


We further have 


N 


N 


rank p^B, + ^ X* = rank (MX*) < rank (M) = 1. (24) 


n=l 


n=l 


N 


N 


Since B is PSD, the matrix (|B-|-^ Pn^n+Yl PnQn] has full rank. Hence, the relation (l24l) 


n=l n=l 
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implies that 


rank fX*) =rank 





X* =1. 


(25) 


n=l 


n=l 


Since X* is rank-one, it can be written as X* = i* (i*)^. The vector i* is thus the solution to (PI). 
This completes the proof. ■ 


Remark 3. As a convex problem, (PI—SDR) can be polynomially solved via an interior-point 
method IfT^ . to arbitrary accuracy. The optimal solution to (PI) is directly obtained from the 
optimal solution to (PI—SDR), without any postprocessing required. 


V. Numerical Results 

As shown in Fig. [2l we consider the setup with N = 5 TX coils and one RX coil, each of which 
has 100 turns and a radius of 0.1 meter. We use cooper wire with radius of 0.1 millimeter for all 
coils. All the coils are horizontal with respect to the xi/-plane. The total resistance of each TX is 
the same as 0.336D. For the RX, its parasitic resistance and load resistance are rp o = 0.336D and 
ri_o = SOD, respectively. The self and mutual inductances are given in Table HI All the capacitors 
are chosen such that the resonance angular frequency is lu = 6.78 x 2tt rad/second [jH. We assume 
that the constraint thresholds 14 = 30\/2 V and An = 5\/2 A. 
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TABLE I: Mutual/Self inductances (/iH) 



TX 1 

TX 2 

TX 3 

TX 4 

TX 5 

TX 1 

5886.8 

0.3565 

0.1253 

0.3565 

0.2984 

TX 2 

0.3565 

5886.8 

0.3565 

0.1253 

0.2984 

TX 3 

0.1253 

0.3565 

5886.8 

0.3565 

0.2984 

TX 4 

0.3565 

0.1253 

0.3565 

5886.8 

0.2984 

TX 5 

0.2984 

0.2984 

0.2984 

0.2984 

5886.8 

RX 

1.6121 

0.00781 

-0.0296 

0.00781 

0.1508 


For comparison, we consider an uneoordinated benehmark seheme of equal eurrent allocation, 
i.e., all TXs carry the same in-phase eurrent whieh ean be adjusted. In partieular, we eompare the 
WPT performanee for the ease with TXs’ eonstraints on the peak voltage and eurrent, and for the 
ease without TXs’ eonstraints, respeetively. We define the effieieney of WPT as the ratio of the 
delivered load power /3o to the total TX power p, i.e., r] = ^. 

Fig. [3] plots the total TX power p and the effieieney p versus the delivered load power (3q. All 
eurves show the feasible and optimal values. For the ease without TXs’ eonstraints, we observe 
that the WPT effieieneies by using magnetie beamforming and by using the benehmark are 87.1% 
and 62.0%, respectively. Hence, the exaet gain of magmatic beamforming is an increase of WPT 
effieieney by 25.1%. 

For the ease with TXs’ eonstraints, we observe that the magnetie beamforming ean deliver more 
power up to 73.6 W to the RX with the effieieney of 74.1%; while the benchmark can deliver 
at most 36.0 W to the RX load with the effieieney of 62.0%. That is, the deliverable power is 
enhaneed by 104.4% by using magnetie beamforming, over the uneoordinated benehmark of equal 
current allocation. Thus, besides the effieieney improvement, another important benefit of magnetic 
beamforming is the enhaneement of the deliverable power, under praetieal TX eonstraints. 

Fig. [3] also shows that the WPT effieieney deereases for 60.5 < /7o < 73.6. To explain the 
decreasing effieieney and obtain insights for magnetic beamforming, we investigate the eases of 
/3o = 60 W and 73.5 W in the following. The optimal eurrents and power of all TXs are obtained 
in Table HIl For = 60 W, TX 1 earries almost the peak eurrent, and most energy is eonsumed 
by TX 1 that has the largest mutual inductanee with the RX. This implies that the TXs with larger 
mutual induetanees with the RX are favorable to earry higher eurrents, to aehieve high effieieney of 
WPT. For = 60 W, all TX eonstraints are inaetive, and it ean be further eheeked that the eurrent 
of each TX is proportional to its mutual inductance with the RX. This verifies (fTSI) in Theorem [T] 
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Si, 
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CD 

o 

& 

X 


■ Mag. B. w.o. TX constr. 


-Benchmark w.o. TX const. 

—B— Mag. B. with TX const. 

— -0^ Benchmark with TX const. 



20 40 60 

RX load power /3o: (W) 


1 


[ ]□□□□□□□□□□□□□□ 


I o.6f<>^<x>0<x>0<>00<)0d^^ 
I 0-4 
^ 0.2 


- Mag. B. w.o. TX constr. 

-- Benchmark w.o. TX const. 

—B— Mag. B. with TX const. 

— -0“ ■ Benchmark with TX const. 


0 20 40 60 

RX load power ( 5 q \ (W) 


80 


100 


Fig. 3: Power and efficiency v.s. RX power /Sq. 


To support higher RX load power of 73.5 W, we observe that TX 1, 3 and 5 carry the (maximally 
allowable) peak current, and TX 2 as well as TX 4 increase their carried currents. The cost is the 
decreased overall efficiency, due to smaller mutual inductance between TX 2, 3, 4, 5 and the RX, 
than that between TX 1 and the RX. 


TABLE II: Optimal results for /3o = 50, 56. 



0 

II 

0 

do = 73.5 


(-7.0698,68.25) 

(-7.071,74.66) 

(* 2 :^ 2 ) 

(-0.0342,0.0016) 

(-3.499,2.22) 

(* 3 :^ 3 ) 

(0.1296,0.023) 

(-7.071,9.62) 

(* 4 ,^ 4 ) 

(-0.0342,0.0016) 

(-3.499,2.22) 

(* 5 :^ 5 ) 

(-0.6617,0.598) 

(-7.071,14.60) 


VI. Conclusion 

This paper has studied the optimal magnetic beamforming for an MRC-WPT system with multiple 
TXs and a single RX. We formulate an optimization problem to minimize the total power drawn 
from the voltage sources of all TXs by jointly designing the currents flowing through different 
TXs, subject to the RX constraint on the required load power and the TXs’ constraints on the peak 
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voltage and current. For the special case of identical TX resistances and neglecting TXs’ constraints, 
the optimal current of each TX is shown to be proportional to the mutual inductance between its 
TX coil and the RX coil. In general, the formulated non-convex QCQP problem is recast as an 
SDP problem. Its SDR is shown to be tight. Numerical results show that magnetic beamforming 
significantly enhances the deliverable power and the WPT efficiency, compared to the uncoordinated 
benchmark scheme of equal current allocation. Furthermore, it is numerically shown that TXs with 
large mutual inductances with the RX are favorable to carry high currents, to achieve high efficiency 
for WPT. 
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